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ABSTRACT 

The search for a luminosity function for galaxies both alternative or compan- 
ion to a Schechter function is a key problem in the reduction of data from catalogs 
of galaxies. Two luminosity functions for galaxies can be built starting from two 
distributions of mass as given by the fragmentation. A first overall distribution 
function is the Kiang function that represents a useful description of the area and 
volume distribution of the Poisson Voronoi diagrams. The second distribution, 
that covers the case of low mass galaxies, is the truncated Pareto distribution 
: in this model we have a natural bound due to the minimum mass/luminosity 
observed and an upper bound ( function of the considered environment) repre- 
sented by the boundary with the observed mass/luminosity overall behaviour. 
The mass distribution is then converted into a luminosity distribution through 
a standard mass-luminosity relationship. The mathematical rules to convert the 
probability density function are used and the two new functions are normalised to 
the total number of galaxies per unit volume. The test of the two new luminosity 
functions for galaxies that cover different ranges in magnitude was made on the 
Sloan Digital Sky Survey (SDSS) in five different bands; the results are compa- 
rable to those of the Schechter function. A new parameter that indicates the 
stellar content is derived. The joint distribution in red-shift and flux , the mean 
red-shift and the number density connected with the first luminosity function for 
galaxies are obtained by analogy with the Schechter function. A new formula 
that allows us to express the mass as a function of the absolute magnitude is 
derived. 



Subject headings: Galaxies: fundamental parameters — Galaxies: statistics - 
Galaxies: luminosity function, mass function 
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1. Introduction 

Over the years the search for a luminosity function for galaxies has played a relevant 
role in the analysis of data from catalogs. A model for the luminosity of galaxies is the 
Schechter function 



$* L 

*(DdL = (-)(-; 



exp( 



L 
~L? 



)dL 



where a sets the slope for low values of L , L* is the characteristi c luminosity and $* is 
the normalisation. This function was suggested by ISchechterl (|1976T ) in o rder to substitute 
other analytical expressions, see for example formula (3) in Kiangl ( 1961 ). Over the years 



this function has also been applied to describe phys ical quantiti e s rela ted to the optical 
luminosity , such as the C O luminos i ty for g alaxies , iKeres et al.l (120031 ). and the barionic 
mass function of galaxies , iBell et al.l (j2003al ). 



An astronomical form of equation (jTJ can be deduced by introducing the distribution 
in absolute magnitude 



$(M)dM = (0.4/nl0)$*10 a4(a+1)(M *- M) 
x exp(-10 a4 ( M *- M ))dM , 



(2) 



where M* is the characteristic magnitude as derived from the data. This distribution has a 
maximum at 



M 



p,max 



M* — 1.085 In (a + 1.0) 



(3) 



where p, max means position of the maximum. In approaching this value the function will 
progressively flatten. 

At present this function is widely used and Table [T] reports , as an example , the 
parameters from three catalogs 



The 2dF Galaxy Reds hift Survey (2dFGRS) based on a sample of 75589 galaxies, see 
first line of Table 3 in iMadewick et al.l J2002h . 



The r*-band luminosity function for a sample of 147 9 86 ga laxies at z = 0.1 from the 
Sloan Digital Sky Survey (SDSS) , see iBlanton et al.l (120031 ). 



The galaxy luminosity functio n for a sampl e of 1 0095 galaxies from the Millennium 
Galaxy Catalogue (MGC) , see briver etaD J2005h . 



Over the years many modifications have been made to the standard Schechter function 
in order to improve its fit: we report three of them. When the fit of the rich clusters 
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lumi nosity function is n o t sati sfactory a two-component Schechter-like function is introduced 



see 



7 Dwarf 



Driver fc Phillippd (119961 ) . This two-component function is the Schechter function when 

* Dwarf dependence when L min < L < L Dwarf : L Dwarf 



< L < L„ 



and has 



J Dwarf 



represents the magnitude where dwarfs first dominate over giants , aowarf the faint slope 
parameter for the dwarf population , the index min and max denote the minimum and the 
maximum. 

Anot her examp l e is th e hybrid Schechter+power-law fit to fit the faint-end of the K- 
band, see bell et all (ho03bh . 



Another function introduced in or der to fit the case of e xtremely low luminosity galaxies 
is the double Schechter function , see iBlanton et al.l (120051 ) , where the parameters $* and 
a that characterise the Schechter function have been doubled in 0* i i,0* i 2 and a.i,a.2- The 
previous efforts bring the attention toward two ranges in luminosity for galaxies : an overall 
zone from high luminosity to low luminosity and the low luminosity zone. This s it uatio n 



remembers the case of t he stars in which th r ee zon es are considered , see IScalol (119861 ) 



Kroupa et al.l (119931 ) , and iBinney fc Merrifieldl (119981 ) ; in this case the range of existence of 
the zones as well as the exponent that characterises the power law behaviour are functions 
of the investigated environment. These three zones in th e mass distr i bution of the stars have 
been investigated at the light of the physical processes in lElmegreenl (120041 ); they correspond 
to brown dwarf masses ~ 0.02.M©, to intermediate mass stars and high mass stars. 

The starting point of this paper is a statistical distribution in the mass of the galaxies 
, A4 as given by a a standard gamma variate with a range of existence < A4 < oo . This 
distribution describes the area of the irregular Voronoi polygons. This distribution in mass 
can be converted in a new statistical distribution for the luminosity of galaxies through an 
analogy with the physics of the stars. This new distribution in luminosity ,L ,is characterised 
by a range of existence < L < oo and a local maximum , named mode. A second 



distribution in the masses starts from a truncated Pareto distribution , after iParetol (118961 ) , 
with a range of existence M. m i n < M. < Ai max . The standard procedure of conversion from 



Table 1: The parameters of the Schechter function from 2dFGRS , SDSS and MGC . 



parameter 
M* [mags] 
a 

$* [h Mpc- 3 } 



2dFGRS 
-19.79 ±0.04 
-1.19 ±0.01 
'1.59±0.1)10- 2 



SDSS (r*) band 
-20.44 ±0.01 
-1.05 ±0.01 

(1.49±0.04)10~ 2 



MGC 
-19.60 ±0.04 
-1.13 ±0.02 
1.77 ± 0.15)10- 2 
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from mass to luminosity allows us to derive a truncated Pareto type luminosity function that 
has a range of existence L min < L < L max . In this distribution the mode is at L min . 

Section 12.21 first introduces the 2D Voronoi diagrams and then describes the mathemat- 
ical details that allow us to deduce two new physical functions for luminosity of galaxies and 
Section [3] reports a first test based on the SDSS photometric catalog. 

In SectionHJthe red-shift dependence of the Schechter function and the first new function 
are explored in detail. Section \5\ reports the mass evaluation for galaxies based on the first 
luminosity function as well as a new formula for the limiting mass. Section [6] summarises 
the results. 



2. From the mass to the magnitude 

The first paragraph briefly introduces the 2D Voronoi diagrams as produced by two 
types of seeds. These two statistical distributions adopted to fit the area of the irregular 
Voronoi polygons can be taken as a starting point to construct two luminosity functions for 
galaxies. 



2.1. The Voronoi Diagrams 



When the seeds are randomly and uncorrelated distributed, are called Poisson Voronoi 
diagrams. A great number of natural phenomena are described by P oisson Voronoi dia- 
we cite some of them : lattices in quantum field theor y . see iDrouffe fc Itzykson 



grams 



(|1984j) ; conductivity and percolation in granular composites , see IJerauld et al.l (ll984al ) and 



Jerauld et al.l (1984 a) ; mod elling growth of metal clusters on amorphous substrates , see 
Dicenzo fc Werthei 5 (11989|); the statistical mechanics of simple glass forming systems in 



2D , see iHentschel et al.l (120071 ); modelling of materi al interface evolution in grain growth 
of polycrystalline materials , see iLee fc Chenl (120061 ). When the seeds are randomly and 



uncorrelated distributed, are called Poisson Voronoi diagrams. A great number of natural 
phenomena are described by Poisson Voronoi diagrams , we cite some of them : lattices 
in quantum field theory , see 



Drouffe fc Itzvksonl ( 



granular composites , see I Jerauld et al. Jl984al ) and 



1984) ; conductivity and percolation in 



Jerauld et al.l (1198401); m odelling growth 



Dicenzo fc Wertheiml 



sec 



Hentschel et al. 



1989); the statistical 



of metal clusters on amorphous substrates , see 
mechanics of simple glass forming systems in 2D 
mater ial interface evolution in grain growth of polycrystalline materials , see ILee &; Chen 



(120071 ): modelling of 



(120061 ) . A review of the Voronoi diagrams applied to the spatial distribution of the galaxies 
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can be found in IZaninettil (120061 ). Here we are interested in the fragmentation of a 2D layer 
of thickness negligible in respect to the main dimension. A typical dimension of the layer 
can be found as follows. The averaged observed diameter of the galaxies is 



D obs w 0.6D 



obs 

max 



Km 

2700 = 27 Mpc 

sec 



(4) 



where D 



obs 
max 



4500 



Km 



corresponds to the extension of the maximum v oid visible on the 



CFA2 slice s. In the framework o f the th eory of the primordial explosions .see lCharlton fc Schramm 
( 119861 ) and IZaninetti fc Ferrarol (119901 ) , this means that the mean observed area of a bubble 
,A obs , is 



j~\obs 

A obs w 47r (^H£)2 = 2290M pc 2 



The averaged area of a face of Voronoi polyhedra , Ay, is 



J^obs 

Ay = 

W 



(5) 



(6) 



sec 



Okabe et al. 


avciagc 

(1992): 


U llUlilUCI 

Zaninetti 


Ul lilt 

(2006) 



dron , L v , is 



v 7 ! 



obis 



12 Mpc 



(7) 



The th ickness of the layer , 6 , can be derived from the shock theory , see iBowers fc Deeming 
(11984J ). and is 1/12 of the radius of the advancing shock , 



S 



D 



obs 



2 x 12 



l.UMpc 



(8) 



The number of galaxies in this typical layer , Ng, can be found by multiplying w 0.1 , 
the density of galaxies , for the volume of the cube of side 12 Mpc : i.e. Ng ~ 172. 

The more common way to insert the seeds of the Voronoi polygons is a random sequence 
in the X and Y direction , see Figure [TJ 

The distribution of the area of the irregular Voronoi polygons is fitted with a Kiang 
function, see formula (1A5I) in Appendix [Al: 



H(x:c) 



Tie) 



[ex 



|C-1 



exp(— ex) 



(9) 



and the captions of Figure [T] also report the value of c as deduced from the parameters of 
the sample. 




Fig. 1. — The Voronoi-diagram in 2D when random seeds are used. The selected region 
comprises 77 seeds and c = 5.1 . 
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Another way to insert the seeds is through a truncated Pareto distribution , see Ap- 
pendix [Bj this is an example of non-Poissonian seeds. 

In polar coordinates the radial distribution of seeds p(r) will vary according to 

P(0°c^ , (10) 
where r is the distance from the center of the box , see Figure [2J 



2.2. A first main physical luminosity distribution 



In order to start, we briefly review how a probability density function ( in the following 
pdf) f(x) changes to g(y) when a new variable y{x) is introduced . We limit ourselves to 
the case in which y(x) is a unique transformation. The rule for transforming a pdf is 



g{y) 



l— I 

I dx I 



We start by assuming that the masses of the galaxies are distributed like a Kiang function, 
see formula ( 1A5I) in Appendix [A] or the previous paragraph. This assumption is justified 
by the fact that the various processes that lead to the formation of a galaxy can follow a 
random fragmentation in 2D. The first transformation is 



x 



M 

M* 



and therefore equation (1A5j) changes to 



V(M)dM 



(M_y 

KM* I 



r c 



-Mr A A 

M* m 

a- 



(12) 



M* 



(13) 



This is a gamma distribution with a scale parameter Ai* and a shape parameter c, and 
its averaged value is 

(M) = cM* . (14) 



The mass- luminosity relationship in the case o f the s t ars is well established both from a 
theoretical point of view, L oc M 3 or L oc M 4 ,see iLangi (119991). a nd fro m an observational 
point of view , L oc 

M 3.43 j n the 

case of MAIN,V , see IZaninettil (120051 ) for further details. 
A power law which is introduced by analogy regulates the relationship between mass and 
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Fig. 2. — The Voronoi-diagram in 2D when truncated Pareto seeds are used. The selected 
region comprises 124 seeds and d = 0.27 as deduced from MLE , see Appendix [B] . 
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luminosity of galaxies, but in this case the regulating parameter a does not have a theoretical 
counterpart. The second transformation is 

M , (15) 



M* \L* 

where 1/a is an exponent that connects the mass to the luminosity. The pdf ( fl3l ) is therefore 
transformed into the following: 

x<kp ("(f)°) <<l • (16) 

where \&* is a normalisation factor which defines the overall density of galaxies , a number per 
cubic Mpc. The mathematical range of existence is < L < oo ; conversely the astronomical 
range is L min < L < L max . 

The relationship connecting the absolute magnitude, M , of a galaxy with its luminosity 

is 

— = ]_o°- 4(Mftoi .©- M) , (17) 



where M bo i tQ is the bolometric luminosity of the sun , which according to ICoxi (120001 ) is 
M fcoZi0 =4.74. 

The third and last transformation connects the luminosity with the absolute magnitude 

^(M)dM = (OAlnlO^— )tf*iO a4( - )(M *- M) 

aT(c) 

xexp(-l0 OA( - M *- M ^)dM . (18) 

This data oriented function contains the parameters M* ,a, c and which can be derived 
from the operation of fitting the observational data. Other interesting quantities are the 
mean luminosity per unit volume, j , 

j = r L9{L)dl = L*vr , (19) 

Jo r ( c ) 



and the averaged luminosity ,(L) 



(L) = f = L* Tic + a) . (20) 
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The density of galaxies is 



J 



and the mean separation between galaxies , 

4 = ™: 1/3 . 



(21) 



(22) 



The symbols j, n* and d* are introduced as in iPadmanabhanl (119961 ). 

Another way to compute the density of galaxies , now n^, of the Ai — L function is 

n*„ = ^* . (23) 

The position of the maximum in magnitudes is at 



M P:max = M* - 1.085 In (c) a . 



(24) 



2.3. The luminosity distribution for low luminosity galaxies 



The Pareto distribution can model nonnegative data with a power law probability tail. 
In many practical app l icatio n s, it is natural to consider an upper bo u nd th at truncates the 
tail lCohen fc Whittenl (jl988h : bevoto k Martnezl Jl998h : Aban et all (bood ). The truncated 
Pa reto distribu t ion h as a wide range of a pplica tions , we list some of them : data analy- 



sis 



Aban et al.l (120061 ) and iRehfeldt et al.l (1l992l ); forest fire area in the Australian Capital 



Territory, fault offsets in the Vernejoul coal field, hyd rocarbon volumes in the Frio Strand 
Plain exploration play and fault lengths on Venus, see Burroughs fc Tebbend (120011 ) . 



In the case of the stars , the low mass distribution of masses , see iSalpeterl (119551 ). can be 
represented by a law of the type p(Ais) oc A4s~ 2 ' 35 , where p(Ais) represents the probability 
of having a mass between A4$ and J\A$ + dM. . By analogy we introduce a truncated Pareto 
distribution, see Appendix [B] , for the mass of galaxies 



* LL (M)dM 



C 



' M W- 
,M* ' 



-A 



M 

M* 



(25) 



where the index LL stands for Low Luminosity and the range of existence is M. m in < -M < 
M-max where min and max denote the minimum and maximum mass. Once the constant C 
is computed as in Appendix [B] we obtain 



(26) 
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Exactly as in the previous case we introduce the transformation represented by equation (fT5l) 
that connects the mass with the luminosity and the distribution in luminosity is 

d-\-a 

9 LL (L)dL = 9l L d Jf\ a , (27) 

\ \ L* ) \ L* I J a 

with the range of existence as L min < L < L max and ^f* LL representing the normalisation. 
The mean luminosity per unit volume, j , 

/ d-\-a d+a \ 

_ t A [ T 2 / Lmax \ a _|_ J 2 / L m j n \ a \ 

j = / L* LL ( L )di = -i L± —^ -V — - ■ ^ 

JL min (^d-a)((^ff-) a - (^ff-) a )L* 

The distribution in magnitude is 

d(M* — am) 

01 (M*-M m „)J a-0 4 ( M *~ M """) d \ 

1(J " — ID a a 



with the range of existence as M m i n < M < M max . This distribution in magnitude contains 
the parameters M min and M max which are the minimum and maximum magnitude of the 
considered catalog and the parameters a , d , and ^ ll which are derived from the fitting of 
the data. 



3. Application to a real sample of galaxies 

The data of the luminosity function for galaxies in five bands of the SDSS are available 



at 
in 



Blanton et all (1200 If ). 



http: / / cosmo.nyu.edu/blanton/lf.html and are discussed from an astronomical point view 



The analysis of the new luminosity function was split in two. The case from high 
luminosities to low luminosities was fitted by \&(M) , equation ffTSl) . The dat a have been 



proces sed through the Levenberg-Marquardt method ( subroutine MRQMIN in iPress et al. 



(Il992h ) in order to find the three parameters a, M* , ty* ; c conversely is introduced by 
hand. In order to associate a statistical probability to each fit we have chosen a range in 
magnitude such as M < M max where M max represents the selected maximum magnitude of 
the sample. 

The results are reported in Table [2] together with the derived quantities j, n* , c?* and 
their uncertainties. Table [2] also reports M max , the number of elements N belonging to the 
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sample , the merit function \ 2 an d the associated p-value that has to b e understood as th e 
maximum probability to obtain a better fitting, see formula (15.2.12) in lPress et al.l (119921 ): 



p = 1 - GAMMQ( 



N-3 
2 ' 2 



(30) 



where GAMMQ is a subroutine for the incomplete gamma function. 
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The uncertainties are found by implementing the error propagation equation (often 
called law of errors of Gauss). The low luminosities range conversely was fitted through 
^ll(M) , equation (T2§j) and the results are reported in Table EJ 



Table 2. Parameters of Fits to Luminosity Function in SDSS Galaxies through the M — L function 



Band 


u* 


g* 


r* 


i* 


z* 


c 


1.1 ± 0.2 


1.0 ± 0.2 


1.1 ± 0.2 


2 ± 0.2 


1.7 ± 0.2 


M* [mags] 


-16.58 ± 0.018 


-18.29 ± 0.008 


-18.77 ± 0.007 


-18.26 ± 0.01 


-18.79 ± 0.004 


** [h Mpc-' A ] 


0.069 ± 0.001 


0.043 ± 0.0003 


0.043 ± 0.000 


0.032 ± 0.0002 


0.034 ± 0.003 


a 


1.40 ± 0.007 


1.32 ± 0.003 


1.5 ± 0.002 


1.74 ± 0.003 


1.70 ± 0.014 


j [mags] 


1.40 L*** 


1.18 L*** 


1.50 L*** 


4.39 L*** 


3.2 L*** 


n* [Mpc -3 ] 


0.097 ± 0.022 


0.051 ± 0.012 


0.066 ± 0.016 


0.14 ± 0.02 


0.11 ± 0.04 


d [Mpc] 


2.17 ± 0.16 


2.67 ± 0.21 


2.47 ± 0.19 


1.91 ± 0.09 


2.06 ± 0.14 


M max [mags] 


-15.78 


-18.2 


-19 


-19.3 


-20 


N 


483 


404 


400 


471 


442 




321 


386 


233 


325 


649 


p = l-GAMMQ(¥f2-,£ 


'-) 


0.31 





1.19 10" 7 


1.0 



I 

I— ' 

I 
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The Schechter function, conversely fits all the range in luminosities and Table H] reports 
the data that come out from the fitting procedure. Table 0] also reports M p ^ max , the value in 
magnitude where the Schechter function peaks ; this value is defined when a — 1 , otherwise 
we leave the box blank. 



Table 3. Parameters of Fits to Luminosity Function in SDSS Low Luminosities Galaxies through the ty LL function 



Band 


g* 


r* 


i* 


z* 


d 


0.3 ± 0.1 


0.4 ± 0.1 


0.5 ± 0.1 


0.9 ± 0.1 


M* [mags] 


-17.2 ± 0.1 


-18.8 ± 0.1 


-17.39 ± 0.1 


-19.3 ± 0.1 


^ LL [h Mpc~ 3 ] 


0.043 ± 0.0043 


0.040 ± 0.0040 


0.026 ± 0.0032 


0.035 ± 0.003 


a 


2.2 ± 0.1 


1.3 ± 0.1 


2.7 ± 0.1 


2.3 ± 0.1 


range [mags] 


-18.2 < M < 16.33 


-19.0 < M < 16.31 


-19.3 < M < 17.22 


-20.0 < M < 17.48 


N 


194 


273 


237 


297 


X 2 


204 


379 


476 


313 



I 

I— ' 

I 
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Table [5] reports the x 2 °f the two zones of the new physical function, their sum and \ 2 
of the Schechter luminosity function. 



Table 4. Parameters of Fits to Luminosity Function in SDSS through the Schechter function 



Band 


u* 


g* 


r* 


i* 


z* 


a 


-0.90 ± 0.01 


-0.88 ± 0.007 


-1.04 ± 0.004 


-0.99 ± 0.005 


-1.07 ± 0.02 


M* [mags] 


-17.92 ± 0.006 


-19.38 ± 0.004 


-20.43 ± 0.003 


-20.81 ± 0.004 


-21.18 ± 0.017 


Mp tm ax [mags] 


-17.92 


-19.38 




-20.81 




$* [h Mpc~ 3 ] 


0.030 ± 0.0003 


0.021 ± 0.0001 


0.015 ± 0.00008 


0.0147 ± 0.00008 


0.0135 ± 0.00006 


j [mags] 


0.95 L*$* 


0.92 L*$* 


1.02 L*<5>* 


0.99 L*<5>* 


1.04 L*$* 


n* [Mpc~ 3 ] 


0.029 ± 0.0003 


0.02 ± 0.0001 


0.015 ± 0.00007 


0.014 ± 0.00008 


0.014 ± 0.00007 


d [Mpc] 


3.24 ± 0.013 


3.64 ± 0.006 


4.02 ± 0.006 


4.08 ± 0.008 


4.12 ± 0.007 


N 


483 


599 


674 


709 


740 


x 2 


330 


753 


2260 


2282 


3245 


p = l-GAMMQ{^,^-) 


5.96 10~ 8 


0.99 


1.0 


1.0 


1.0 
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The Schechter function , the new function and the data are reported in Figure [3] , 
Figure H] , Figure [5] , Figure [6] , and Figure [7] when the u*,g* ,r* , i* and z* bands are 
considered; Figure [S] , Figure M , Figure [TD] , Figure [TT] , and Figure [£2] report the residuals 
of the u*,g* ,r* , i* and z* band. We have used H = lOO/i km s _1 Mpc^ 1 , with h — 1 in 
all the numerical evaluations. Due to the testing phase of the new M. — L function, we have 
o mitted the propagatio n of other values of h on the derived quantities , see the discussion 
Blanton et all J200lh . 



m 



Table 5. Synoptic \ 2 ■ 



Band 


u* 


9* 


r* 


V 


z* 


X 2 physical luminosity function 


321 


386 


233 


325 


649 


X 2 luminosity function low luminosities 





204 


379 


476 


313 


X 2 sum of two zones 


321 


590 


612 


801 


962 


X 2 Schcchtcr luminosity function 


330 


753 


2260 


2282 


3245 
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» I i i i I i i i i I i i i i I i i i i I i i i i I i I 

-20 -10 -18 -T7 -16 

absolute magnitude M 

Fig. 3. — The luminosity function data of SDSS(u*) are represented through the error bar. 
The fitting continuous line represents our two luminosity functions ( ffl8l) and (129]) ) and the 
dotted line represents the Schechter function. 
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» L i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

1 -22 -21 -20 -19 -18 -17 

absolute magnitude M 

Fig. 4. — The luminosity function data of SDSS(g*) are represented through the error bar. 
The fitting continuous line represents our two luminosity functions ( (fIBj) and (129]) ) and the 
dotted line represents the Schechter function. 
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_l 



■9 I i .... i .... i .... i .... t .... i .... t ... . 

1 -23 -22 -21 -20 -19 -18 -17 

absolute magnitude M 

Fig. 5. — The luminosity function data of SDSS(r*) are represented through the error bar. 
The fitting continuous line represents our two luminosity functions ( (fIBj) and (129]) ) and the 
dotted line represents the Schechter function. 
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» i .. i .... i .... i .... i .... i .... i ... . 
-23 -22 -21 -20 -19 -16 

absolute magnitude M 

Fig. 6. — The luminosity function data of SDSS(i*) are represented through the error bar. 
The fitting continuous line represents our two luminosity functions ( (fIBj) and fl29|) ) and the 
dotted line represents the Schechter function. 
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Fig. 7. — The luminosity function data of SDSS(,2*) are represented through the error bar. 
The fitting continuous line represents our two luminosity functions ( (fIBj) and (129]) ) and the 
dotted line represents the Schechter function. 
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Fig. 8. — The residuals of the fits to SDSS(w*) data. The empty stars represent our two 
luminosity functions ( (fl8|) and (|29|) ) and the filled points represent the Schechter function. 
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Fig. 9. — The residuals of the fits to SDSS(g*) data. The empty stars represent our two 
luminosity functions ( ([18]) and (|29]) ) and the filled points represent the Schechter function. 
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Fig. 10. — The residuals of the fits to SDSS(r*) data. The empty stars represent our two 
luminosity functions ( (fl8|) and (|29|) ) and the filled points represent the Schechter function. 
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Fig. 11. — The residuals of the fits to SDSS(i*) data. The empty stars represent our two 
luminosity functions ( ([18]) and (|29|) ) and the filled points represent the Schechter function. 
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Fig. 12. — The residuals of the fits to SDSS(z*) data. The empty stars represent our two 
luminosity functions ( (ITS]) and (12U|) ) and the filled points represent the Schechter function. 



31 



The value obtained for the parameter a should be compare d with that of the normal 
stars which is 3 or 4 as suggested by the theory, see for example iLana (119991) . or w 3.8 as 
suggested by the observations for M. > 0.2A^ Q , see for example I Cox! (120001 ) . When the 
three classes of stars are considered we have a =3.43 (MAIN), a =2.79 (GIANTS) and a=2.43 



(SUPERGIANTS), see IZaninettil ( I2005f ) for further details. 



The variation of j when the range in magnitude is finite rather than infinite can be 
evaluated by coupling together formula (TTTT) and ([TBI 



-15.78 



1O O.4(A/ 6oi , -M)^ M ) dM 



(31) 



21 



On inserting the parameters of SDSS band u* ( which is the case in which the M. — L 
function covers all the range in magnitude of the data , see Tabled) and M 6o ; j0 =M u . =6.39 
, j = 1.4 10 8 Lq is obtained. This value increases by 5.63 % when the ran ge is infini t e; see a 
similar discussion concerning Schechter's function around formula (33) in lLin et al.l (119961 ). 



In absence of observational data that represent the luminosity function, we can generate 
them through Schechter's parameters, see Table [TJ Thi s is done, for example for the 2dF 



Galaxy Redshift Survey (2dFGRS) .see lCross et al. 



(120011 ) and data of the Schechter function 



in Tabled! The parameters of the A4—L function are reported in Table[6]where the requested 
errors on the values of luminosity are the same as the considered value. 



4. Tests involving z 

Some useful formulae connected with the Schechter function in a Euclidean ,non-relativistic 
and homogeneous universe are reviewed ; by analogy new formulae for the first part of the 
M. — L function are derived. 



Table 6: The parameters of the M. — L luminosity function 
based on 2dFGRS data ( triplets generated by the author) 

2dFGRS 

c 0.1 

M*[mags] -19 ±0.1 

^*[h Mpc~ 3 } 0.4 ±0.01 

a 1.3 ±0.1 
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4.1. The behaviour of the Schechter function 



The flux of radiation , / , is introduced 



./ 



47rr 2 



(32) 



where r represents the distance of the g alaxy . The joint distribution in z and / for galaxies 
, see formula (1.104) in iPadmanabhanl (119961 ) , is 



dN 



dfldzdf 



4tt( 



^crit 



(33) 



where dVt , dz and df represent the differential of the solid angle , the red-shift and the 
fl ux respectively. The L* of difference between the previous formula and formula (1.104) 
in 



Padmanabhanl (119961 ) is due to the small difference in the definition of $. 
The formula for z cr u is 



^crit 



4^7c| 



(34) 



where cl represents the light velocity; the CODATA recommen ds cj, = 299792.458^ . The 



mean red-shift of galaxies with a flux / , see formula (1.105) in IPadmanabhanl (119961 ) , is 

T(3 + a) 



^crit 



r(5/2 + a) 



(35) 



he n umber density of galaxies per unit flux interval, see formula (1.106) in iPadmanabhan 
(119961 ) , is 



dN 
dlnf 



$* L* . 
T ( 4^7* 



3/2 



r(- + «) 



(36) 



The number of galaxies in z and / as given by formula (1331) has a maximum at z = z. 
where 

Z-max Zcrit^J & ~\~ 2 

The value of z max can be derived from the histogram of the observed number of galaxies 
expressed as a function of z . For practical purposes we analysed the 2dFGRS data release 



(37) 



available at the web site: http://msowww.anu.edu.au/2dFGRS/. In particular we added 



together the file parent.ngp.txt that contains 145652 entries for NGP strip sources and the 
file parent.sgp.txt that contains 204490 entries for SGP strip sources. Once the heliocentric 
red-shift was selected we processed 219107 galaxies with 0.001 < z < 0.25 . A comparison 
between the observed and theoretical number of galaxies as a function of z is reported in 
Figure Q3J 
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Fig. 13. — The galaxies of the 2dF Galaxy Redshift Survey with 15.08 < bJmag < 15.81 
( with bJmag representing the relative magnitude used in object selection), are isolated in 
order to represent a chosen value of / and then organised in frequencies versus heliocentric 
redshift , (empty circles); the error bar is given by the square root of the frequency. The 
theoretical curve generated by the Schechter function of luminosity (formula (1331) and param- 
eters as in column 2dFGRS of Tabled]) is drawn (full line). The theoretical curve generated 
by the M. — L function for luminosity ( formula ( 1381) and parameters as in column 2dFGRS 
of Table [6]) is drawn (dashed line); \ 2= 320 for the Schechter function and x 2= 283 for the 
Ai — L function. 
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Another interesting catalog is the 6dF Galaxy Survey that has measured ar ound 150000 



redshi fts and 15000 peculiar velocities from galaxies over the southern sky , see iJones et al. 



(120061 ). It is available at the following address http: / /vizier. u-strasbg.fr/viz-bin/VizieR?-source=VII/249 
; we selected the re-calibrated bJ magnitude and the recession velocity cz. Figure 



| reports 

the observed and theoretical number of galaxies as a function of z for the 6dF Galaxy Survey. 
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Fig. 14. — The galaxies of the 6dF Galaxy Survey with 14.15 < bJmag < 14.9 ( with 
bJmag representing the relative magnitude used in object selection), are isolated in order 
to represent a chosen value of / and then organised in frequencies versus redshift , (empty 
circles); the error bar is given by the square root of the frequency. The theoretical curve 
generated by the Schechter function of luminosity (formula (|33l) and parameters as in column 
2dFGRS of Table [T]) is drawn (full line). The theoretical curve generated by the M. — L 
function for luminosity ( formula (1381) and parameters in column 2dFGRS of Table [6]) is 
drawn (dashed line); x 2= 1373 for the Schechter function and x 2= 1197 for the M. — L 
function. 
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4.2. The behaviour of the M. — L function 

The joint distribution in z and / , in presence of the M. — L luminosity (equation ffT6lT 



is 



dN c , 5 4 . z 2 



47r(— )V*(^-) . (38) 

The mean red-shift is 



2 g+c . , 2 c+3 a 

, Z)=Z J_ 4-^-r ( 2a+ C 2^ 

s 7 T(c + 3/2a) v ; 



The number density of galaxies per unit flux interval is 

dN 1 L* 3/2 ^T (c + fa) 

~ 16 7r 3 / 2 / 3/2 r (c) ' 

The number of galaxies as given by formula ( |38j) has a maximum at -z ma x where 



J crit 



(40) 



; C + a) a/2 . (41) 



A comparison between the observed and theoretical number of galaxies as given by the M. — L 
function is reported in Figure [13] where the 2dF Galaxy Redshift Survey is considered and 
in Figure [H] where the 6dF Galaxy Survey is considered. 



5. Mass evaluation 

One method to deduce the mass of a star by its absolute visual magnitude is presented; 
the mass of a galaxy is deduced by analogy. In the case of the galaxies, the bolometric 
correction of the stars will be replaced by the sun's absolute magnitude and mass-luminosity 
ratio different in each selected band. 



5.1. The case of the stars 



In the case of the stars it is possible to parameter i se the mass of the star , Ais , as a 
function of the observable colour (B — V) , see lZaninettil (120051 ) . The first equation connects 
the (B — V) colour with the temperature 

2~bv 



(B — V) — K BV + 



T 



(42) 



here T is the temperature, K^y and Tbv are two parameters that can be derived by im- 
plementing the least square method on a series of calibrated data. The second equation 
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describes the bolometric correction , BC, 

BC = M hol - My = -10 log 10 T + K BC , (43) 

where M ho \ is the absolute bolometric magnitude, My is the absolute visual magnitude, T B c 
and are two parameters that can be derived through the general linear least square 
method applied to a series of calibrated data. The third equation is the usual formula for 
the luminosity 

logio(^4 = 0-4(4.74 -M bol ) , (44) 

where L is the luminosity of the star and L Q the luminosity of the sun. The fourth equation 
is the usual mass-luminosity relationship for stars 

L M s 
logioll— ) = °lm + b hM log 10 (— — ) (45) 

L,q JVl Q 

for M > 0.2M Q , 
where M.s is the mass of the star and Ai® is the mass of the sun. 
With these four equations the mass of the star is 

. Ms 

-0.4 My - 0.4 K BC + 4.0 In ( (b _%1 Kbv ) (in (lO)" 1 ) 



4 T BC (-(B-V)+K BV ) + Lgg6 _ ^ 



(46) 



with the various coefficients as given by Table 1 in IZaninettil (120051 ). As an example, the 
mass of a star belonging to MAIN SEQUENCE V is 

■° s '°£r- 7 - 769+0 - 8972 ln ( (fi -J) 6 + 1 o.64ii ) • (47) 

We can now express the colour (B — V) as a function of the absolute visual magnitude 
My and the following formula for the mass of the star is obtained 

( 9378 \ 

l°g,o j± = -7.769 + 0.8972 In ( - {9me _^ mMv) ) (48) 

MAIN SEQUENCE V when - 5.8 < My < 11.8 , 
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where W is the Lambert W- function, after iLambertl (117581 ). A test of the previous formula 

can be done at the two boundaries: when My =-0.58 , log 10 = 1.63 against the calibrated 

value log 10 j^ 2 - = 1.6 and when M y =11.8 , log 10 xrr = —0.56 against the calibrated value 

lno- ^S- 
1O Sl0 M a 



—0.66, see Table 3.1 in lBowers fc Deeming! (I1984I ). 



5.2. The case of the galaxies 

The mass of a galaxy can be evaluated once the mass luminosity ratio , R is given 



R 



<f> 



(49) 



, R = 20 bv 


Padmanabhan 


( 


1996) 


Bell & de Jong 


(2001 


), demonstrated 



R = 5.93 by 



(119611) and iPersic fc Salucci Jl992h 
van der Marell (119911 ). Further on 



M 



colour, and therefore, type. If the bright end of the luminosity function is dominated by mas- 
sive, evolved, red galaxies, and the faint end by low mass, blue galaxies, then 4^ oc L~ 0M 
(GIANTS III) at the bright end and ^ oc L ~°- 7 (MAIN SEQUENCE V) at the faint end , 



see coefficients of Table 1 in IZaninettil (120051 ). Then ^ will almost certainly not be constant 



due to different prevailing po pulations of star s at th e boundaries of the luminosity function. 
The scatter in the models by iBell fc de Jong! (120011 ) is a starting point when evaluating the 
validity of assuming a constant Generally, near-infrared ^ ratios are more constant than 
optical passband, but still vary with luminosity. In our framework , we made R a function 
of the passband , in order to have the same results for the masses of the galaxies once the 
absolute magnitude is given, see Figure [151 In our framework R can be expressed as 



R 



(M) 
(L) 



(50) 



On inserting formula (fT4"j) and formula (TSOj) in the previous ratio the following formula for 
M* is found 

T(c + a)M Q 



M* = RL* 

r( c ) l 

From equations (1151) and (1511) a formula for the mass of the galaxy is found 



(51) 



M 



BMP 



AM, 



bol,Q 



-0.4M*) r(c + fl) ( 10 -0. 



4M+0.4M*\ a 



cT (c) 



(52) 



An application of the previous formula is reported in Figure ( 1151) , where the mass of 
galaxies as a function of the absolute magnitude in the five bands of SDSS is drawn. In 
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this Figure Mhni n is different for ea ch selected band and equal to the value suggested in 
equation (16) of iBlanton et al.l (120011 ). 
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Fig. 15. — Logarithm of the mass of the galaxy as a function of the absolute magnitude. The 
SDSS bands are u* with M feoij0 =M M * =6.39 and R = 6 (full line) , g* with M bol:Q =M g * Q =5.07 
and R = 13 (dashed), r* with M boZj0 =M r » =4.62 and R — 16 (dot-dash-dot-dash), i* with 
M bo i :Q =Mi* Q =A.52 and R = 15 (dotted), z* with M 6o , i0 =M 2 » =4.48 and i? = 14 (dash-dot- 
dot- dot). 
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The new formula fl52l) allows us to deduce the mass of the galaxy from its absolute 
magnitude and can be easily particularized in different pass-bands. As an exa mple , with the 
data o f SDSS in the five bands reported in Tabled, M bo ^ Q as in equation (16) of lBlanton et al. 
(120011 ) and R as in Figure [T51 we have 



M= 215234 e- 0S579M M Q 
M= 97216 e-°- 6978M A4 

M = 490000 e-° m41M M & 

M= 2691000 e-°- 5294M A4 

M = 3434000 e-°- 5175M M Q 



u* band when - 20.6 < M < -15.7 
g* band when - 22.0 < M < -18.2 
r* band when - 23.0 < M < -19 
i* band when - 23.5 < M < -19.3 
z* band when - 23.7 < M < -20.0 



(53) 



The method here suggested to deduce the mass of the galaxy can be compar e d wit h the 



formula that com es out from the Tully-Fisher relation , see iTully fc Fisherl (119771 ) and 
Tully et all (I1998I ). In the Tully-Fisher framework the mass of a rotating galaxy can be 
parameterised as 



M = 50V f 4 A4 e 



Km 4 

Km 



(54) 



sec 



McGaughl ( 120051 ). The mass to 



where Vf is the rotational velocity expressed in 
light ratio in our framework scales oc L^ 1 ^^" 1 with a depending on the selected catalog and 
band. This ratio oscillates , referring to the SDSS data, between a minimum dependence 
in the g* band , ^ oc L~ 24 and a maximum d ependence in the i* b and , ^ oc L~ 0A2 . 



A comparison should be made with ^ oc L 35 in Ivan der Marell (119911 ) for a sample of 37 



bright elliptical galaxies ; this result was obtained by implementing axisymmetric dynamical 
models. The completeness of the mass sample of the galaxies belonging to a given catalog 
can be evaluated in the following way. The limiting apparent magnitude is known, and is 
differe nt for each catalog. In the case of the SDSS (r* band) is m=17.6, see iBlanton et al. 

feooih . 



The corresponding absolute limiting magnitude is computed and inserted in equa- 
tion (!52j) . The limiting mass for galaxies , M. L ,is 



lo (O.4M 6o! , -O.4M*) r ( c + a ) ( 10 -0Am+2.0Logl0(^) +10.0+0.4 M* 



Mi 



R- 



c(T(c)) 



-M 



5 



(55) 



where clz is the radial distance expressed in km/s. In order to see how the parameter z 
influences the limiting mass, Table [7| reports the range of observable masses as a function of 
z. 



6. Conclusions 



We have split the analysis of the luminosity function in two . The analysis of the main 
new luminosity function , formula (|18p . from low luminosities up to maximum magnitude 
shows that , see Table |2l 



1. The parameter c varies between 0.1 and 2. It must be remembered that the theory 
predicts 2 , 4 and 6 for the ID, 2D and 3D fragmentation respectively. 

2. Parameter a varies between 1.32 and 1.74. The numerical mass-luminosity relationship 
for the stars gives values of the parameter a comprised between 2.43 and 3.43. 

3. The M. — L function represents a better fit of the observational data in comparison 
with the Schechter function once the concept of maximum magnitude of the sample is 
introduced. Without this limiting magnitude the situation is inverted. 



The case of low luminosities galaxies was describe by a truncated Pareto type luminosity 
function , see formula ( |29l) . This new luminosity function is described by two physical 
parameters , d , and a denoting respectively the distribution in mass and the mass-luminosity 
connection. The analysis of the data for low luminosities galaxies as reported in Table [3] shows 
that 



1. The parameter d varies between 0.3 an d 0.9 . This ya lue should be compared with d 
of the stars which is 2.3 — 1 = 1.3 , see iKroupal (1200 ll ). 



2. The parameter a varies between 1.3 and 2.7. 



Table 7: The limiting mass for the SDSS catalog ,tt* band, when R and Mt> [ :& are those of 
Figure [TBI The limiting apparent magnitude is m=17.6. 



mass range z 

1A110 8 M Q <M<1.7 10 U M Q 0.001 

3.791O 9 A^ <M<1.7 10 U M Q 0.01 
lM0 n M e <M<1.7 10 u M e 0.1 

1.410 n .M„ <M<1.7 IQ^Mv 0.13 
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The theoretical number of galaxies as a function of the red-shift presents a maximum 
that is a function of a and / for the Schechter function and c , a and / for the first M. — L 
function ; the agreement with the maximum in the observed number of galaxies is acceptable. 

The observable range in masses can be parameterised as a function of z and the ratio 
between maximum and minimum luminosity is 232 at z =0.001 but drops to 1.1 at z =0.15 , 
see Table [0 

Perhaps a more comprehensive way of comparing the mass estimates of equation (j54j) 
( Tully-Fisher relation ) with those given here (equation fl52|) and equation (1531) ) would 
be as follows. We take the same sample of galaxies for which the luminosity function was 
computed in Section [3] and compute their mass function according to a given value of R . 
This mass function can be compared with those of other galaxies in a common passband and, 
also in this case, the distribution is expressed through a Schechter function. Three cases are 
now analysed 



1. The Bell case , see Bell et all (12003af ). where $* = O.OlMpcr 3 /Log w {M), M* 



5.3 ltt m M G) and a = -1.21. 



2. The Bottema case , see iBottemal (119971 ). where $* = O.OUMpc' 3 / Log 10 (M), M* 



2.24 l^Mr, and a = -1.20. 



3. Kennicutt-Kroupa case , see iKennicuttl (119831 ); iKroupa et al.l (119931 ) . where $* 
0.01lMpc- 3 /Log w (M), M* = 3.78 10 10 M Q and a = -1.22. 



Figure [161 reports the already cited standard distributions as well as our \& {M.) / 'Logio(M) 
when the range in masses is that given by the conversion from luminosity to masses. 



-44- 




Fig. 16. — Barionic Mass function of galaxies: our ty(Ai)/ Logio(Ai) (full line) , Bell case 
(dashed), Bottema case (dot-dash-dot-dash) and Kennicutt-Kroupa case (dotted). In our 
case we considered the SDSS band u* with Mfe o ; i0 =M tt * Q =6.39 and R=6.0 . 
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The analysis of the two new functions for the luminosity for galaxies here derived gives 
a marginally better fit but certainly the Schechter function for it's simplicity and fewer 
parameters still represents a good model for the luminosity function for galaxies. At present 
the study of the Schechter function is not yet terminated and two new equations were derived: 
equation ([3]) that represents the maximum in magnitude distribution and equation (|33|) that 
gives the value of z at which the observed number of galaxies is maximum. 



A. On the Kiang Function 

The starting point is the distribution in length , s , of a segment in a random fragmen- 
tation 

p(s) = A exp (— \s)ds , (Al) 

where A is the hazard rate of the exponential distribution. Given the fact that the sum , u 
, of two exponential distributions is 

p(u) = A 2 wexp (— Xu)du . (A2) 

The distribution of ID Voronoi segments , I, ( the midpoint of the sum of two segments) can 
be found from the previous formula by inserting u = 21 

p(l) = 2Xlexp(-2Xl)d(2Xl) . (A3) 

On transforming in normalised units x = v we obtain 

p(x) = 2xexp(-2x)d(2x) . (A4) 



Afte r one century o f studies on the Voronoi diagrams , see the two memories IVoronoil (119071 ) 



and IVoronoil (119081 ) , the law of the segments in ID is the unique analytical result on the 
field. When this result is expressed as a gamma variate we obtain formula(5) of lKian gl (jl~96fih 



H(x;c) = (cxY -1 exp(— cx) , (A5) 
r(c) 

where < x < oo , c >0 and T(c) is the gamma function with argument c; in the case 
of ID Voronoi diagrams c = 2. It was conjectured that the area in 2D and the volumes 
in 3D of the Voronoi diagrams may be approximated as the sum of two and three gamma 
variate of argument 2. Due to the fact that the sum of n independent gamma variates with 
shape parameter q is a gamma variate with shape parameter c = Y^t c ii the area and the 
volumes are supposed to follow a gamma variate of argument 4 and 6. This hypothesis 
was later named "Kiang's conjecture", and the equation ( 1A5I) used as a fitting function , 
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see Kumar et al.l (119921 ); IZaninettil (120061 ). or as a n hypothes is to a c cept or to reject using 
the standard procedures of the data analysis, see iTanemural (119881 . 120031 ). A new way to 
parametrise the ID, 2D and 3D cells on the base of th e considered dimensionality has been 
introduced , see formula (12) in lFerenc &: Nedal (120071 ). 



B. On the Truncated Pareto Distribution 



The sta r ting p df (probability density function) is the Pareto distribution iParetol (118961 ); 
Evans et all (boom . P, 



P(x; a, c) 



ca 



x 



c+1 



where a < x < oo , a >0,c > 0. The average value is 

_ ca 

x 



c- 1 

which is defined for c > 1, and the variance is 



a 2 c 



a 



(c-2)(c-l 



,2 ' 



(Bl) 



(B2) 



(B3) 



which is defined for c > 2. The presence of an upper bound ,b, allows us to introduce the 
following pdf , named truncated Pareto Pt , 



Pt(x; a, b, c) 



1 ca c 
1 - (f) c x c+1 



(B4) 



here a < x < b , a > ,b > 0, b > a and c > 0. The distribution function of the 
truncated Pareto is 



F(x; a, b, c) 



l"(f) c 



The average value of the truncated Pareto pdf is 



x 



ca i-qr 1 

c-l l-(f) c 



and the variance of the truncated Pareto pdf is 

„ numerator 



a 



with 



denominator ' 
denominator = (c - 2) (c - l) 2 (a 2c - 2 b c a c + b 2c ) 



(B5) 



(B6) 



(B7) 



:bs) 
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numerator = 

cb 2 a 2c + 2 a 2+c c 2 b c - 4 cV +1 6 c+1 + 2 c 2 b 2+c a c - a 2+c cb c 
+2 c 3 a c+1 6 c+1 - a 2+c c 3 6 c - c 3 b 2+c a c - cb 2+c a c + ca 2 b 2c . 

The variance of the truncated Pareto is always defined for every value of c > 0; conversely 
the variance of the Pareto distribution can be defined only when c > 2. The parameter c 
can be derived through the maximum likelihood estimator (MLE). The likelihood function 
is defined as the probability we would have obtained a given set of observations if given a 
particular set of values of the distribution parameters, q, 

L( c ) = f(x! . . .x n \ci . . .c n ) . (B9) 

If we assume that the n random variables are independently and identically distributed, then 
we may write the likelihood function as 

n 

L(c) = f(x 1 \c 1 ...c p )...f(x n \c 1 ...c p ) = Y[f(x i \c 1 ...c p ) . (BIO) 

i=i 

The maximum likelihood estimates for the q are obtained by maximising the likelihood 
function, L(c). Equivalently, we may find it easier to maximise lnf(xi) , termed the log- 
likelihood. So, for a random sample X\ . . . x n from a truncated Pareto distribution, the 
likelihood function is given by 

n 

L( C ) = Uc^r 1 - (&TT 1 (^ c+1 )" 1 • (Bn) 

i=l 

In this model we have assumed that a— min(xi . . . x n ) and b= max(xi . . . x n ). 
Using logarithms, we obtain the log-likelihood 

n 

lnL(c) = ncln(a) + nln -r— — In Xj . (B12) 

1 ~ {v° i=i 

Taking the first derivative , we get 

|>L(c) = 

nlna + !! + !!iil!Ml_^ ln ,, = o . (bis) 

C 1 ^b) i= i 

The parameter c can be found by solving numerically the previous non-linear equation. 



I thank Massimo Ramella , Tao Kiang and the anonymous referee for comments. 
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